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Abstract 

Vibrational  energy  transfer  between  an  IR  laser,  a  solid  surface  and 
an  adatom  is  studied  by  a  generalized  master  equation  which  includes  memory 
effects.  Numerical  problems  associated  with  the  temporally  delocalized 
memory  kernel  are  overcome  by  introducing  the  isomnesic  approximation. 

Results  obtained  for  time  scales  ranging  from  2  ps  to  2  us  for  H(H2P)/KCi(001 ) 
show  that  effective  vibrational  excitation  is  possible  with  low  power  lasers 

p 

(1  W/cm  )in  spite  of  fast  phonon  relaxation.  Detailed  time-evolution  of  the 
system  is  seen  to  differ  radically  from  that  predicted  by  Markovian  theories, 
even  for  fairly  long  time  scales.  - 


I.  Introduction 

Developments  in  the  experimental  study  of  clean,  well -characterized 
solid  surfaces  and  surface  processes  such  as  scattering,  adsorption,  de¬ 
sorption  and  migration  in  the  presence  of  laser  radiation  have  resulted 
in  important  technological . applications  (e.g.,  laser-assisted  vapor 
deposition  of  microcircuits).  The  underlying  dynamics  of  these  hetero¬ 
geneous  systems  is  a  many-body  problem  involving  the  elementary  excitations 
of  the  solid  and  its  surface,  such  as  phonons,  electrons  and  plasmons,  the 

degrees  of  freedom  of  the  adsorbed  species  or  gas  molecules,  and  those  of 
1  2 

the  laser.  ’  The  excitation  process  can  be  considered  as  direct  laser¬ 
pumping  of  the  substrate  leading  to  a  primarily  non-selective  thermal 
effect,  or  a  pumping  of  the  adsorbed  species,  which  can  be  selective  in 
the  sense  that  only  a  single  degree  of  freedom  is  excited  by  the  laser. 

The  resulting  energy  absorption  and  the  dissipation  involving  all  other 

3  8 

degrees  of  freedom  govern  the  overall  chemical  process. 

In  this  work,  the  vibrational  excitation  of  a  prototype  system,  viz. 
the  adsorptive  bond  (adbond)  between  an  atom  and  a  surface  is  studied. 

This  system  is  not  in  the  same  category  as  those  in  which  internal  vi¬ 
brational  modes  of  an  adsorbed  molecule  are  excited,  with  subsequent  energy 

7  9 

transfer  to  the  adbond.  *  It  is  substantially  simpler,  in  view  of  the 
aim  of  this  work,  which  is  to  investigate  the  mechanisms  involved  in  IR 

laser-stimulated  surface  processes  (LSSP)  using  a  treatment  which  is  as  close 
to  first-principles  as  possible.  One  of  the  fundamental  questions  addressed 
is  whether  a  low-power  laser,  because  of  its  coherence,  can  deposit  large 
amounts  of  energy  into  a  specific  mode  of  a  system  (such  as  an  adbond) 
which  is  coupled  to  energy  sinks,  such  as  phonons,  which  are  not  infinite. 
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Early  classical  treatments^  of  the  problem  and  phenomenological  quantum 
4  5 

treatments  ’  provided  conflicting  answers  to  this  question.  Recent,  more 

g 

fundamental  studies  have  suggested  that  extremely  high  laser  powers  would 
be  needed  to  effectively  desorb  an  adspecies.  However,  our  preliminary 
work^  showed  that  a  very  small  degree  of  detuning  between  the  laser  and 
energy  levels  of  the  system  excluded  the  possibility  of  long-term  energy 
absorption.  The  reason  for  this  behaviour  is  seen  to  be  the  oscillatory 
nature  of  the  laser  energy  absorption  "rate"  (actually,  the  memory  kernel) 
for  the  adbond,  with  a  frequency  related  to  the  detuning,  typically  a  small 
percentage  of  the  laser  frequency  itself,  viz.  -  lO^-lO11  s"1.  Thus,  for 
processes  lasting  over  a  few  hundred  picoseconds,  substantial  cancellations 
in  the  overall  absorption  can  be  expected,  and  only  a  very  high-power  laser 
would  be  able  to  compete  with  the  loss  of  energy  from  the  adbond  to  phonons 
through  relaxation  effects.  In  view  of  this,  the  systems  chosen  here  have 
little  or  no  detuning,  viz.  systems  with  one  bound-to-bound  transition 
exactly  resonant  with  the  laser  frequency  and  another  very  close  to  or 
equal  to  the  Debye  frequency  of  the  solid.  This  represents  an  ideal 
situation  which  emphasizes  the  coherent  nature  of  laser  radiation. 

The  widely  varying  temporal  behaviour  of  dynamical  systems  far  from 
equilibrium  is  well  established.  From  a  fundamental  point  of  view, 
theoretical  treatments  of  such  systems  depend  crucially  on  the  time  scales 
being  considered.  The  most  common  long-time  limit,  for  example,  allows 
a  description  of  energy  transfer  in  terms  of  average  relaxation  rates  in 
the  Markovian  approximation,  where  approach  to  equilibrium  is  assumed  to 
be  monotonic.  Non-Markovian  behaviour  is  assumed  to  last  for  extremely 
short  times  compared  to  typical  relaxation  times  of  the  system.  The  forces 
are  considered  to  be  random,  with  a  very  short  correlation  time.  Resonant 
excitation  with  a  low-power  laser  does  not  satisfy  these  conditions,  and 


one  should  expect  non-Markovian  behaviour  to  persist  for  times  that 

increase  with  decrease  in  laser  power.  Whether  the  two  approaches  lead  to 

the  same  steady  state  is  also  open  to  question. 

Non-Markovian  effects  are  conveniently  treated  by  using  a  generalized 

master  equation  (GME)  to  describe  the  time-evolution  of  the  adbond  (the 

"relevant"  system)  with  the  laser  and  phonon  fields  as  two  uncoupled 

irrelevant  systems  (reservoirs).12"1^  The  Liouville  equation  for  such  a 

system  is  conveniently  solved  by  using  the  Zwanzig  projection-operator 
12 

technique  and  results  in  a  set  of  coupled  integrodifferential  equations 

14 

for  the  occupation  probabilities  P^  of  states  |S>  of  the  adbond: 

PSU)  =  l  J*  dt*  [Kss,(t-f)PSi(t')  -  Ks,s(t-t')Ps(t')]  ,  (1) 

S‘*S  0 

where  the  memory  kernels  Kss,  include  contributions  due  to  the  phonon  field 
and  the  radiation  field  K^: 

Kss,(t)  =  Kss*)(t)  +  Kss')(t)-  (2) 

Using  a  one-dimensional  lattice  and  a  Debye  model  to  describe  the 
lattice  vibrations,  we  have  obtained  in  a  closed  form  which  involves 

a  summation  over  all  lattice  positions  £  and  number  of  phonons  n  inde¬ 
pendently.1^’^  The  former  is  due  to  a  summation  of  pair  potentials 
between  the  adatom  and  the  lattice  atoms,  while  the  latter  results  from 
an  expansion  of  exp[yu, (t)]  in  a  Taylor  series,  where  y££l(t)  is  the 
lattice  displacement  correlation  function 

yu.(t)  =  <<u£(t)u£.»  .  (3) 

Using  techniques  described  elsewhere,11  y££i(t)  can  be  obtained  in  closed 
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where  R  ,  and  I  ,  Involve  products  of  polynomials  In  t  and  cos  (wQt) 
or  sin  (w pt ) ,  and  Is  the  Debye  frequency.  The  phonon  part  of  the 
memory  kernel  can  then  be  written  In  terms  of  all  the  system  parameters  as 
follows: 


^(t,‘Ap WI&  <w‘)n 

X  X£^uu'^S'^EnHi,*^t^C0S(wcc't)  +  ^  sin(uw,t)]. 


,D  (°)  2 

AP  a  8  frV> 


Bp  =  6  (0(O)kbT)2/ms-h(uo3  , 


E"“'(t)  =  j  .0^4...  R“'6(t)  ’ 

6 

<J#-l)/2  n-i  j 

W(t)  =  j  4,3,5..  V0!(n-j0)!  RU^t}  TU,{t)  ’ 


max(je,j0)  <  n 


8^  Is  a  smallness  parameter  which  dictates,  in  part,  the  convergence 

properties  of  the  expansion  in  number  of  phonons  n.  It  is  of  the  order 

of  .02  for  common  temperatures  T  and  common  Debye  frequencies  (kg  is  the 

Boltzmann  constant  and  m$  the  mass  of  a  lattice  atom),  y  and  y'  take  on 

the  values  1  and  2.  The  supermatrix  X  involves  powers  of  exp(a)  and 
o 

exp(c<  u  >>)  and  the  matrix  elements 

.  .  u(z-zj0b 

E(“!  .  <S|e  |S'>  < 


where  a  is  the  lattice  constant,  u  is  the  lattice  displacement  operator, 

and  2  is  the  distance  of  the  adatom  from  the  equilibrium  position  of  the 

outermost  lattice  atom.  The  frequencies  10^,  correspond  to  transitions 

betweea  levels  |S>  and  | S ' >  of  the  adbond  represented  by  a  Morse  potential 

with  depth  ,  equilibrium  distance  and  exponent 

( r ) 

The  radiation  part  of  the  memory  kernel  is  given  by 

( r) 

K$$.(t)  *  Arzs,s  zs$,  cos  |)t, 

2  2 

where  ZSS '  =  <slzlsl>*  The  amplitude  =  (IeQ  /egCh  )  is  proportional 
to  the  laser  frequency  I  and  the  square  of  the  charge  separation  e^  between 
the  adspecies  and  the  surface. 

The  complicated  form  of  R££l  and  I  ,  precludes  analytic  solution  of 
the  6ME.  Numerical  solution  is  made  difficult  by  the  presence  of  high 
frequencies  in  and  the  large  difference  in  the  magnitudes  of 

Ar(-1017  s'2  for  a  1  W/cm2  laser)  and  Ap(-1029  s“2  for  =  0.1  eV). 

One  can  only  solve  the  exact  GME  for  a  few  periods  of  oscillation  of  K^,  , 
which  amounts  to  a  few  picoseconds. 

Examination  of  a  typical  phonon  kernel.  Fig.  1,  shows  that  it  has 
substantial  amplitude  only  for  times  of  the  order  of  a  Debye  period  x^  , 
representative  of  phonon  correlation  times.  The  laser  kernel,  on  the 
other  hand,  has  constant  amplitude  for  the  case  ,  for  a  given 

pair  of  levels  ( S ,S ' ) .  For  times  larger  than  xp,  the  integrated  effect 
due  to  will  be  unchanging  because  of  cancellations,  whereas  that 

due  to  K<.^  will  be  cumulative.  We  have  recently  demonstrated  this  effect^ 
using  the  exact  GME  (with  a  very  high  laser  power  in  order  for  the  effect 
to  manifest  itself  in  the  very  short  time  for  which  solution  was  possible), 
and  predicted,  via  extrapolation  arguments,  that  effective  pumping  of  the 
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adbond  with  a  low-power  laser  would  be  possible  If  the  long-time  regime 
could  be  examined.  We  present  here  an  actual  calculation  of  this  long¬ 
time  behaviour  using  low-power  lasers,  and  find  that  the  average  adbond 
energy,  given  by 

E(t)  =  |  Es  P$(t)  (13) 

where  are  eigenvalues  of  the  zeroth  order  adbond  Hamiltonian,  rises 
rapidly  after  an  initial  loss  to  the  phonon  field,  and  exhibits  a  decidedly 
non-Markovian  pattern  throughout. 

The  approximations  required  for  the  long-time  treatment  are  described 
in  Section  II.  The  resulting  probability  and  energy  profiles  are  presented 
in  Section  III,  where  comparison  is  made  with  the  exact  results  (for  a  short- 
time  scale)  and  with  a  Markovian  limit. 


II .  The  Isomnesic  (Constant-Memory)  Approximation 

We  begin  with  the  assertion  that  the  physically  important  situation  is 
one  where  the  laser  frequency  and  the  Debye  frequency  equal  different 
transition  frequencies  the  adbond.^  For  this  case  K^|(t)  becomes 

a  constant  kss,  (from  the  term  with  [w^i  I  55  p^us  0SCll1ating  terms  which 
we  ignore  in  the  sense  of  the  rotating-wave  approximation.^  The  phonon 
term,  Kj|?^(t),  has  a  complicated  behaviour  typified  by  Kg^(t)  in  Fig.  1, 
but  we  note  that  it  has  appreciable  magnitude  for  times  =  2  ps,  and  almost 
vanishes  subsequently.  If  the  mesh  size  of  our  theoretical  experiment  is  no 
less  than  2  ps,  we  can  treat  Ks^(t)  as  a  delta  function,  so  that,  from 
Eq.  (2), 

K$$(  (t)  =  Ps$l  6 (t)  +  k^1  .  -  (14) 


'SS 


where 


flee,  =  /  dt  k4?}(t) 


(15) 
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In  Eq.  (14),  the  first  term,  originating  from  the  coupling  of  the  adbond 
to  the  phonons,  is  the  amnesic  or  Markovian  term,  and  would  lead,  in  the 
absence  of  the  laser  term,  to  the  elimination  of  all  memory  effects.  The 
second  term,  due  to  the  resonant  laser  coupling,  is  the  isomnesic  non- 
Markovian  term,  and  represents  a  constant  memory  effect  all  the  way  back 
to  t=0. 

The  subsequent  treatment  is  simplified  by  defining  "diagonal"  elements 
nSS  as 


n 


SS 


.1  * 


S'*S 


SS' 


(16) 


and  similarly  for  k^ 
and  matrix  notation 
written  as 


,  and  by  introducing  the  vector  notation  JP(t)  =  { Ps ( t ) } 
a  =  {figgil  »  k  =  {k^,}  .  etc.  The  exact  GME  may  be 


P(t)  = 


0 


t 

/  dt '  K(t-t')  •  P(t')  ; 


(17) 


in  the  isomnesic  approximation  (IA)  it  becomes 

t 

P(t)  =  k  •  /  dt 1  P(t')  +  n-P(t)  . 

=  0 


(18) 


Formally  this  can  be  solved  without  any  further  approximations.  Differentiating 
both  sides  of  Eq.  (18)  we  get  the  second-order  differential  equation 

P(t)  =  k-P(t)  +  n  -P(t)  ,  (19) 


whose  Laplace  transform, 

M(s).  P(s)  =  I(s)  ,  (20) 

where 

M(s)  e  k  +  sn  -  s2a  ,  (21) 

T(s)  E  -  sP(0)  .  (22) 

and  P(s)  is  the  Laplace  transform  of  ^(t),  is  easily  inverted  using  the 
Heaviside  expansion  theorem. 

While  this  effectively  eliminates  the  non-Markovi an  bottleneck,  it  is 
instructive  to  generate  solutions  of  a  Markovian  version  of  this  class  of 
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problems.  To  do  this,  we  have  to  make  provisions  for  localizing  KSS '  1n 
time.  This  brings  to  the  fore  the  pervasive  issue  of  time  scales  in 
dynamical  many-body  systems,  our  specific  question  being  the  relative 
degree  of  localization  ascribable  to  and  .  If  the  laser 

radiation  is  localized  as  a  pulse,  it  can  still  extend  over  a  period  of  time 
much  longer  than  the  extent  of  K^|^(t).  In  other  words,  the  delta  function 
representing  K<-^(t)  may  be  quite  different  from  the  delta  function 
representing  K^,  (t),  and  the  former  may  in  fact  not  really  be  considered 
as  localized  with  respect  to  the  latter.  With  these  provisos  in  mind,  we 
can  proceed  to  write  the  Markovian  version  of  Eqs.  (14),  (17),  (19)  and  (20): 


K(t)  «  (n  +  f)  <5 ( t )  ,  (23) 

P(t)  =  W  •  P(t)  ,  (24) 

J(s)  •  P(s)  =  U  ,  (25) 

where 

f  =  k  t  ,  (26) 

W  =  n  +  f  ,  (27) 

J(s)  =  (n  +  f  -  sa)  ,  (28) 

U  =  -  P(0)  ,  (29) 


and  t  represents  the  actual  extent  of  the  laser  signal  or  a  variable  time 
parameter.  Equation  (24)  is  the  well-known  Pauli  master  equation  (PME). 

Ill .  Results  of  the  Exact,  Isomnesic  and  Markovian  Approximations 
The  formalism  of  Sections  I  and  II  was  applied  to  the  system 
H( H^O )/KC2.(001 ) ,  the  primary  reasons  for  choosing  which  are  the  shallow 
adsorption  potential  well  and  the  small  number  of  bound  states  (7)  it 
generates.17  In  Fig.  2  the  results  of  the  IA  and  the  exact  GME  are  compared  for 
the  system  of  references  11  and  15.  The  two  are  identical  for  H(H20)/KU,  and 
are  not  shown.  While  this  is  a  short-time-scale  comparison,  the  closeness 


of  the  two  leads  us  to  expect  the  long-time  IA  results  to  be  reliable. 

The  Markovian  approximation,  on  the  other  hand,  cannot  provide  an  un¬ 
ambiguous  comparison  because  of  large  variations  in  the  results  for 
different  values  of  t  ,  shown  in  Fig.  3.  All  subsequent  figures,  there¬ 
fore  present  our  IA  results. 

Figure  4  is  a  composite  display  of  the  function  E(t)  obtained  from 

results  of  Eq.  (28)  for  six  ranges  of  time  differing  by  factors  of  10  and 

2 

for  a  laser  intensity  I  =  1  W/cm  .  It  is  seen  that  in  going  from  a  pico¬ 
second  time  range  to  a  micro-second  time  range,  the  behaviour  of  the  system 
changes  dramatically.  The  early  behaviour  is  dominated  by  fast  energy 
transfer  to  phonon  modes,  completely  overwhelming  the  slow  laser  pumping. 

The  monotonic  decay  is  representative  of  relaxation  phenomena.  For  the  100  ps 
range,  a  false  steady  state  seems  to  have  been  reached,  but  is  in  fact  only 
an  artifact  of  the  large  difference  between  Ar  and  Ap.  The  adbond  energy 
starts  to  increase  in  the  ns  range,  again  suggesting  a  possible  steady  state. 
However,  this  is  a  transitional  period  during  which  the  relative  importance 
of  the  phonon  terms  and  the  laser  terms  starts  to  change.  The  latter,  in 

Ct  t 

fact,  originate  from  a  probability  profile  of  the  form  (Ae  cos  gt  - 
Beat  sin  gt),  where  a  +  ig  are  complex  roots  of 

det  M(s)  =  0  ,  (30) 

and  A  +  iB  are  complex  coefficients  obtained  from  M(s)  and  T(s). 

The  phonon  terms  come  from  real  roots  and  are  therefore  nonosci 1 latory. 
The  magnitudes  of  the  exponents  a  in  the  laser  part  of  the  solution  are 
much  smaller  than  corresponding  exponents  in  the  phonon  term.  This  is 
responsible  for  the  radically  different  nature  of  E(t)  in  the  different 
time  regimes.  Thus,  the  seemingly  monotonic  rise  in  E(t)  apparent  in 
the  10  ns  regime  is  seen  to  be  only  the  early  segment  of  a  cycle  which 
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Figure 

Fig.  1 
Fig.  2 
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Fig.  4 
Fig.  5 


Captions 


.  Memory  kernels  K^l(t)  for  (S,S‘)  *  (0,3).  Solid  line:  3-*0;  dotted 
line:  (K3. 

.  Probability  profiles  for  system  of  references  11  and  15  obtained  by 
solving  the  exact  6ME  (solid  line)  and  the  isomnesic  6ME  (dashed  line). 

.  Average  adbond  energy  E( t)  obtained  within  the  Markovian  approximation 
for  a  range  of  values  of  the  parameter  t. 

.  E(t)  obtained  by  using  the  isomnesic  approximation  for  time  scales 

2 

ranging  from  picoseconds  to  microseconds,  for  laser  intensity  I  =  1  W/cm  . 

2 

.  Same  as  Fig.  4  for  I  =  25  W/cm  . 
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result,  the  sharp  transitions  assumed  here  will  be  broadened  by  a  detuning- 
type  effect,  and  estimates  of  the  desorption  rate  presented  above  will 
have  to  be  modified.  However,  the  selectivity  represented  by  the  dominant 
Ti  energy-transfer  process  will  be  retained. 

For  a  complete  study  of  desorption,  transitions  to  the  continuum 
must  be  included,  and  at  the  outset  it  is  difficult  to  predict  which  levels 
will  be  most  important.  Work  in  this  direction  is  in  progress,  and 
represents  a  departure  from  the  first-principles  approach  in  the  inclusion 

of  the  bound- to-continuum  transitions.  This,  as  well  as  a  complete  first- 
principles  treatment  of  vibrational  LSSP  including  desorption,  will  be 
presented  in  the  near  future. 
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becomes  evident  in  the  1G0  ns  results  and  beyond.  The  frequency  of  this 
oscillation  is  seen  to  be  kg  ,  where  kg  is  the  amplitude  of  the  laser 
kernel  k^,  ,  and  is  essentially  the  Rabi  frequency  for  the  pair  of  levels 
resonant  with  the  laser. 

While  the  frequency  of  the  long-term  oscillations  of  P(t)  and  E(t) 
depends  only  on  the  laser  term,  being  proportional  to  the  square  root  of 
the  intensity,  the  rate  of  onset  of  this  behaviour  depends  on  o  which  varies 

linearly  as  kg  and  inversely  as  the  amplitude  flg  of  the  phonon  kernel. 

2  2 

Thus,  increasing  the  laser  power  from  1  W/cm  to  25  W/cm  advances  the  onset 
time  for  the  oscillatory  state  while  increasing  the  frequency  by  a  factor 
of  5,  as  displayed  in  Fig.  5. 

Physically,  the  results  of  Figs.  4  and  5  can  be  interpreted  in  terms 
of  desorption  rates.  Thus,  the  time  required  for  ] A [ e  and  |B|e  to 
become  large  enough  so  that  the  amplitude  of  E(t)  in  the  oscillatory  regime 
is  substantially  larger  than  the  (t  =  0)  equilibrium  Boltzmann  value, is  a 
reasonable  measure  of  the  desorption  time  to  within  an  order  of  magnitude. 
This  time  is  estimated  by  Tg  =  a”^ .  For  our  case, 

V1  *  a  =  .  5  =  (104  w)  s'1  ,  (31) 

2  2 

where  w  is  laser  power  in  W/cm  .  Thus,  for  a  1  W/cm  laser,  a  hydrogen 
atom  will  stay  adsorbed  on  the  KC 2. ( H ^0 )  surface  for  -  10-4  s  on  the 
average  before  desorbing,  or  at  least  becoming  very  highly  excited.  The 
data  mentioned  here  is  for  T  =  150  K. 

It  must  be  pointed  out  that  a  number  of  mechanisms  whereby  the  desorp¬ 
tion  cross  section  could  be  modified  are  not  included.  Spontaneous  decay, 
phase  relaxation,  phonon-phonon  interaction  and  electron-hole  pair  creation 
are  some  of  the  ones  being  cons)  :red  in  extensions  of  this  work.  As  a 


DL/413/83/01 

056/413-2 


ABSTRACTS  DISTRIBUTION  LIST,  056/625/629 


Dr.  J.  E.  Jensen 
Hughes  Research  Laboratory 
3011  Malibu  Canyon  Road 
Malibu,  California  90265 

Or.  J.  H.  Weaver 

Department  of  Chemical  Engineering 
and  Materials  Science 
University  of  Minnesota 
Minneapolis,  Minnesota  55455 


Or.  W.  Knauer 

Hughes  Research  Laboratory 
3011  Malibu  Canyon  Road 
Malibu,  California  90265 

Dr.  C.  B.  Harris 
Department  of  Chemistry 
University  of  California 
Berkeley,  California  94720 


END 

filmed 

6-85 

DTIC 


